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A  SIMPLE  TECHNIQUE  FOR  DISPLAYING  WDB-II  DIGITAL 
CARTOGRAPHIC  IMAGERY  IN  KEE 


1  Introduction 

The  Knowledge  Engineering  Environment  (KEE)  is  a  sophisticated  commercial 
software  development  environment  by  IntelliCorp,  Inc.,  Mountain  View,  California.  KEE 
contains  a  computer  graphics  toolkit  called  KEEpictures  for  building  graphical  software  to 
be  used  with  KEE.  KEEpictures  is  an  object-oriented  computer  graphics  application  that  is 
fully  integrated  with  the  KEE  system  and  provides  the  developer  with  a  wide  range  of 
capability  from  rendering  simple  displays  to  implementing  complicated  animation 
techniques  [1]. 

Software  systems  developed  using  KEE  often  have  a  requirement  to  display  maps. 
It  is  not  hard  to  integrate  cartographic  databases  with  KEE,  and  there  is  an  abundance  of 
commercial  databases  to  choose  from.  The  Central  Intelligence  Agency  (CIA)  and  Defense 
Mapping  Agency  (DMA)  are  two  excellent  sources  of  digital  map  data  within  the  defense 
community,  for  example.  This  paper  will  describe  a  simple  technique  for  displaying  maps 
drawn  from  the  World  Data  Bank  II  cartographic  database  obtained  from  the  CIA.  The 
maps  were  displayed  on  a  Sun  workstation  in  KEE  using  KEEpictures. 


2  Description  of  the  Hardware  and  Software  Environment 

A  Sun  3/260  workstation  running  Sun  operating  system  (Sun  OS)  version  4.1  was 
used  as  the  development  platform  for  this  work.  KEE  version  3.1.1 17  was  used  with  Sun 
Common  Lisp  version  2.1.3,  and  this  version  of  KEEpictures  supported  only  monochrome 
graphical  operations.  The  draw  program  described  in  section  6  was  written  using  the 
SunCore  computer  graphics  software  library.  SunCore  was  a  software  product  by  Sun 
Microsystems  that  had  been  bundled  with  the  Sun  OS  prior  to  version  4.0.  Sun  has 
discontinued  support  for  the  SunCore  product  in  Sun  OS  version  4.0  and  later  releases, 
and  consequently  SunCore  libraries  and  include  files  are  no  longer  distributed.  These  files 
were  ported  from  the  Sun  OS  3.5  environment  to  the  4.1  environment  without  any 
noticeably  deleterious  side  effects. 


3  Approach 

The  technical  approach  is  described  briefly  in  the  following  five  steps,  and  figure  1 
shows  tnc  processing  pipeline. 

1 .  Read  the  geographic  areas  of  interest  from  files  op  World  Data  Bank  II  (WDB-II) 
distribution  tapes  and  convert  the  files  from  WDB-II  file  format  to  tiles  in  a  simpler 
file  format,  called  sff  files. 

2 .  Optionally  filter  or  correct  the  data. 

3.  Draw  the  map  on  tne  Sun  console. 

4.  Save  a  screendump  of  the  image  in  a  file. 
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5. 


Convert  the  raster  image  into  a  lisp  expression  that,  when  loaded  using  the  KEE 
lisp  listener,  creates  a  Common  Windows  bitmap  of  the  image.  Write  this  lisp 
expression  to  a  file. 

6.  Load  the  file  in  KEE  using  the  lisp  listener,  creating  a  Common  Windows  bitmap 
structure.  The  map  can  be  viewed  in  a  KEEpictures  viewport  by  making  the  bitmap 
structure  the  value  of  the  bitmap  slot  in  a  KEEpictures  bitmap  unit. 


4  World  Data  Bank  II 

World  Data  Bank  II  is  a  cartographic  database  of  the  world  distributed  by  the 
Central  Intelligence  Agency  (CIA),  Washington,  DC.  WDB-D  is  a  digital  representation  of 
the  world  divided  into  five  geographic  areas.  Each  area  is  contained  on  a  single  magnetic 
tape  volume  and  each  tape  contains  three  files  of  unique  line  segments  defined  by 
individually  digitized  points.  The  five  geographic  areas  are  North  America,  South 
America,  Europe,  Africa,  and  Asia,  and  the  three  files  contained  in  each  tape  volume  are  a) 
coastlines,  islands,  and  lakes,  b)  rivers,  and  c)  international  boundaries.  In  addition,  the 
North  America  volume  contains  a  fourth  file  of  internal  boundaries  to  further  delineate  the 
individual  states  of  the  USA  and  the  provinces  of  Canada.  WDB-II  was  digitized  at  scales 
ranging  from  1:1  million  to  1:4  million;  the  nominal  scale  is  1:3  million.  WDB-II  contains 
approximately  6  million  digitized  points.  [2] 


4.1  Description  of  WDB-II  Files 

Each  of  the  16  WDB-II  files  consist  of  a  sequence  of  line  segments  and  are 
formatted  on  the  tape  in  the  following  way.  Each  line  segment  begins  with  a  20  character 
record: 
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Columns 


1  -  7  Line  segment  number 

8  -  9  Rank 

10-15  Number  of  points  in  line 

16-20  Zeroes 

Each  point  in  the  line  segment  follows  and  they  are  described  by  this  20  character 

record: 


Columns  1-7 
8-15 
16-20 


Latitude  (DDMMSSI) 
Longitude  (DDMMSSI) 
Point  number 


Each  line  segment  in  WDB-II  has  a  numeric  rank  assigned  to  it  which  is  described 

below: 


Coast,  islands,  and  lakes 

01  Coasts,  islands,  and  lakes  that  appear  on  all  maps 

02  Additional  major  islands  and  lakes 

03  Intermediate  islands  and  lakes 

04  Minor  islands  and  lakes 

06  Intermittent  major  lakes 

07  Intermittent  minor  lakes 

08  Reefs 

09  Major  salt  pans 

10  Minor  salt  pans 

13  Major  ice  shelves 

1 4  Minor  ice  shelves 

15  Glaciers 

Rivers 

0 1  Permanent  major  rivers 

02  Additional  major  rivers 

03  Additional  rivers 

04  Minor  rivers 

05  Double-lined  rivers 

06  Major  intermittent  rivers 

07  Additional  intermittent  rivers 

08  Minor  intermittent  rivers 

10  Major  canals 

1 1  Canals  of  lesser  importance 

12  Irrigation  type  canals 

International  boundaries  or  limits  of  sovereignty 
0 1  Demarcated  or  delimi  ted 

02  Indefinite  or  in  dispute 

03  Other  line  of  separation  or  sovereignty  on  land 

Internal  boundaries 

0 1  First  order  administrative 
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4.2  Reading  WDB-II  Files  from  Tape  and  Converting  them  to  sff  Files 

The  WDB-II  data  files  tape  format  specification  was  not  very  convenient  for  our 
purposes  so  a  special  program  was  used,  called  wdbtosff,  to  convert  the  data  into  a  form 
that  is  easier  to  use.  The  wdbtosff  program  reads  a  WDB-II  file  and  converts  it  into  an 
ASCII  text  file  called  an  sff  file  (which  stands  for  simpler  file  format)  given  by  the 
following  grammar: 


<sff  file> 

<segment  list> 
<segment> 

<point  list> 
<point> 

<rank> 

<n  points> 
<separator> 
ccharacter  integer> 
<digits> 

<digit> 

<lat> 

<lon> 

ccharacter  float> 


<segment  Iist>  EOF 
<segment>  [<segment  list>] 

<rank>  < separator  <n  points>  <separator>  <point  list> 
<point>  [<point  list>] 

<lat>  <separator>  <lon>  <separator> 

<character  integer> 

<character  integex> 

Space  I  CarriageRetum  [<separator>] 

[+ 1  -)  <digits> 

<digit>  [<digits>] 

0111213141516171819 
ccharacter  float> 
ccharacter  float> 
ccharacter  integer>  [.  cdigits>] 


The  value  of  cn  points>  describes  the  number  of  instances  of  cpoint>  in  cpoint 
list>  which  follow. 

WDB-II  files  can  be  read  directly  from  the  distribution  tape  using  the  UNIX  dd 
command  and  piped  to  the  wdbtosff  program,  whose  output  is  redirected  to  a  file.  The 
following  examples  show  this  in  detail  for  each  WDB-II  file.  Note  that  files  on  the  Asia 
volume  have  a  block  size  of  8000  bytes;  all  other  files  have  a  block  size  of  3600  bytes. 
Note  also  that  WDB-H  tapes  are  written  in  the  EBCDIC  character  set  and  must  be  converted 
to  ASCII  using  the  conv=ascii  option. 


Vol  1:  North  America 

dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  na.cil 
dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  na.riv 
dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  na.bdy 
dd  if=/dev/nxmt0  bs=3600  conv=ascii  I  wdbtosff  >  na.pby 

Vol  2:  South  America 

dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  sa.cil 
dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  sa.riv 
dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  sa.bdy 

Vol  3:  Europe 

dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  eur.cil 
dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  eur.riv 
dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  eur.bdy 

Vol  4:  Africa 

dd  if=/dev/nrmt0  bs=3600  conv=ascii  1  wdbtosff  >  af.cil 
dd  if=/dev/nrmt0  bs=3600  conv=ascii  I  wdbtosff  >  af.riv 
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dd  if=/dev/nrmtO  bs=3600  conv=ascii  I  wdbtosff  >  af.bdy 
Vol5:  Asia 

dd  if=/dev/nrmtO  bs=8000  conv=ascii  I  wdbtosff  >  as.cil 
dd  if=/dev/nrmtO  bs=8000  conv=ascii  I  wdbtosff  >  as.riv 
dd  if=/dev/nrmtO  bs=8000  conv=ascii  I  wdbtosff  >  as.bdy 

The  source  code  of  the  wdbtosff  program  is  shown  in  Appendix  A. 


5  Optionally  Filtering  the  Data  or  Correcting  Boundary  Errors 

A  data  reduction  program  called  reduce  is  shown  in  Appendix  B  and  can  be  used  to 
filter  data  points  against  three  criteria:  minimum  and  maximum  distance  between  two 
points  and  angle  described  by  three  points.  Roughly,  an  adjacent  point  is  discarded  if  the 
distance  between  the  test  point  and  adjacent  point  does  not  lie  within  a  range  specified  by 
the  user,  or,  if  the  angle  described  by  three  sequential  points  is  near  1 80  degrees  then  the 
middle  point  is  discarded.  A  precise  description  of  the  algorithm  can  be  found  in  the  source 
code  listing  in  Appendix  B. 

WDB-n  data  point  values  aie  bounded  by  [90S,  90N]  in  latitude  and  [180W,  180E] 
in  longitude.  The  files  as.riv,  as.cil  and  sa.cil  contain  some  line  segments  that  cross  the 
international  date  line,  which  :nay  cause  problems  for  unsophisticated  drawing  programs 
(such  as  the  one  described  next).  If  a  point  that  crosses  the  date  line  is  not  accounted  for 
correctly  a  line  connecting  a  pair  of  points  on  either  side  of  the  date  line  will  be  drawn 
erroneously  across  the  map.  A  simple  program  called  bcor  corrects  this  problem  by 
splitting  the  offending  segment  into  new  segments,  each  of  which  reaches  the  date  line  but 
does  not  cross  it  The  source  code  listing  for  bcor  is  given  in  Appendix  C. 


6  Drawing  the  Map  and  Making  a  Screendump  of  the  Image 

Maps  are  drawn  on  the  Sun  console  with  a  program  called  draw,  which  projects  the 
data  points  using  a  quasi-equirectangular  projection.  Here  is  the  usage: 

draw  mIMIc  rankmax  latmin  Ionmin  latmax  lonmax  filel  [file2  file3  ...] 

Draw  takes  several  arguments,  the  first  of  which  is  a  single  character  that  identifies 
the  type  of  display  console  to  the  program.  The  three  choices  are  m  which  stands  for  low 
resolution  (1152x900  pixel)  monochrome,  M  for  high  resolution  (1600x1280  pixel) 
monochrome,  or  c  for  low  resolution  (1152x900  pixel)  color.  Segments  with  rank  value 
greater  than  rankmax  are  not  drawn.  The  next  four  arguments  define  two  points  on  the 
globe  which  describe  an  area  that  is  projected  onto  the  display  console:  latmin,  Ionmin, 
latmax,  and  lonmax.  Latitudes  south  of  the  equator  and  longitudes  west  of  the  prime 
meridian  are  represented  as  negative  floating  point  numbers;  all  others  are  represented  as 
positive  floating  point  numbers.  The  arguments  that  follow  are  names  of  sff  files  that  are  to 
be  drawn  on  the  console.  For  example,  to  draw  the  whole  world  on  a  high-resclution 
monochrome  console  the  following  command  line  would  be  given: 

draw  M  15  -90  -180  90  180  na.cil  na.riv  na.bdy  na.pby  sa.cil  sa.riv  sa.bdy  eur.cil 
eur.riv  eur.bdy  af.cil  af  riv  af.bdy  as.cil  as.riv  as.bdy 
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To  draw  European  coastlines,  islands,  lakes,  rivers,  and  international  boundaries  of 
rank  1-3  projected  from  25N  30W  to  75N  35E  on  a  low  resolution  color  console  the 
following  command  line  would  be  given: 

draw  c  3  25  -30  75  35  eur.cil  eur.riv  eur.bdy 

When  the  drawing  is  finished  the  image  should  be  saved  in  a  file  using  the 
screendump  command: 

screendump  >  theworld.ras 

Note:  the  draw  program  should  not  be  run  from  the  console  because  feedback  from 
draw  will  appear  on  the  screen  and  ruin  the  image.  The  source  code  listing  for  draw  is 
shown  in  Appendix  D. 

Figures  2-5  show  several  example  drawings.  Figure  2  shows  all  WDB-II  data. 
The  draw  command  used  here  was: 

draw  M  15  -90  -180  90  180  box2  box3  af.bdy  af.cil  af.riv  as.bdy  as.cil  as.riv 
eur.bdy  eur.cil  eur.riv  na.bdy  na.cil  na.pby  na.riv  sa.bdy  sa.cil  sa.riv 

Box2  and  box3  are  small  Files  that  contain  data  for  drawing  rectangles.  Box2  was 
used  to  draw  the  rectangle  enclosing  the  whole  map.  Note  the  small  rectangle  around  the 
Persian  Gulf  region.  Box3  was  used  to  draw  this  rectangle.  Figure  3  shows  an  enlarged 
view  of  this  region.  The  contents  of  box2  and  box3  are  shown  in  Appendix  F. 

Figure  3  shows  the  Persian  Gulf  region  bounded  by  23N  42E  and  32N  60E.  The 
draw  command  used  to  make  this  picture  was: 

draw  M  15  23  42  32  60  box3  box4  af.bdy  af.cil  af.riv  as.bdy  as.cil  as.riv 

Box3  was  used  to  draw  the  rectangle  enclosing  the  map,  and  box4  was  used  to 
draw  the  small  rectangle  in  the  Kuwait  region.  Figure  4  shows  an  enlarged  view  of  this 
region.  The  contents  of  box3  and  box4  are  shown  in  Appendix  F. 

Figure  4  shows  the  Persian  Gulf  region  bounded  by  23N  42E  and  32N  60E.  The 
draw  command  used  to  make  this  picture  was: 

draw  M  15  23  42  32  60  box4  af.bdy  af.cil  af.riv  as.bdy  as.cil  as.riv 

Box4  w?s  used  to  draw  the  rectangle  enclosing  the  map.  The  contents  of  are 
shown  in  Appendix  F. 

Figure  5  shows  a  less  cluttered  view  of  coastlines,  islands,  lakes,  and  international 
boundaries  of  the  world,  and  internal  boundaries  of  North  America.  The  draw  command 
used  to  make  this  picture  was: 

draw  M  1  -90-180  90  180  box2  af.bdy  af.cil  as.bdy  as.cil  eur.bdy  eur.cil  na.bdy 
na.cil  na.pby  sa.bdy  sa.cil 

Box2  was  used  to  draw  the  rectangle  around  the  map.  The  contents  of  box2  are 
shown  in  Appendix  F. 
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Figure  2. 
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Figure  3. 
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Figure  5. 
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Making  a  Common  Windows  Bitmap  File 


A  program  called  suntokee  converts  a  raster  image  saved  as  a  monochrome 
rasterfile  into  a  file  that  contains  a  lisp  expression  which,  when  loaded  by  the  KEE  lisp 
listener,  assigns  a  Common  Windows  bitmap  structure  containing  the  identical  raster  image 
to  a  lisp  variable.  The  variable  name  bitmappname  is  specified  as  an  argument  to  the 
suntokee  program: 

suntokee  bitmappname  <  rasterfile  >  lispfile 

To  define  a  Common  Windows  bitmap  named  kuwait-bitmap  in  the  file  called 
kuwait-bitmap.lisp  from  the  screendump  image  in  the  1600x1280  rasterfile  kuwait.ras,  for 
example,  the  following  line  would  be  used: 

suntokee  kuwait-bitmap  <  kuwait.ras  >  kuwait-bitmap.lisp 

The  lisp  source  file  is  a  simple  ASCII  text  file  and  can  be  viewed  using  an  editor  or 
printed  out  by  the  curious  user.  Also,  the  file  can  be  compiled  using  the  lisp  function 
compile-file  for  greater  efficiency. 

The  suntokee  program  is  shown  in  Appendix  E. 


8  Viewing  the  Map  in  KEE 

The  next  step  is  to  load  the  lisp  file  with  the  KEE  lisp  listener,  using  the  lisp  load 
function.  Having  created  a  file  called  kuwait-bitmap.lisp  using  the  example  in  section  7, 
the  following  line  would  be  entered  in  the  KEE  lisp  listener 

(load  "kuwait-bitmap.lisp") 

This  creates  a  1600x1280  Common  Windows  bitmap  structure  named  kuwait- 
bitmap.  To  view  the  bitmap  first  make  a  KEEpictures  viewport  programmatically  with 
(unitmsg  'viewports  'create.instance!  ...),  or  interactively  using  the  mouse  and  menu. 
Next  create  a  KEEpictures  bitmap  unit  in  the  viewport.  The  bitmap  unit  contains  a  slot 
called  bitmap ,  which  always  contains  a  reference  to  a  Common  Windows  bitmap  structure. 
Set  the  value  of  the  bitmap  slot  of  this  KEEpictures  bitmap  unit  to  kuwait-bitmap 
programmatically  with  (put. value ...),  or  interactively  using  the  mouse  and  menu.  Figure  6 
shows  the  raster  image  from  the  rasterfile  in  the  KEEpictures  viewport. 

Since  the  geometry  of  the  cartographic  projection  and  the  geographic  boundaries  of 
the  bitmap  image  are  known  it  is  a  straightforward  job  to  implement  coordinate  system 
transforms  from  geographic  coordinate  space  to  KEEpicture  coordinate  space  and  back. 
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Appendix  A 


wdbtosff  Program  Source  Code  Listing 

/'  WDBTOSFF:  a  simplistic  program  that  reads  a  file  in  WDB-II  format 
from  stdin  and  writes  ASCII  text  in  the  following  format  to  stdout: 

<segmentl> 

<segment2> 


where  each  segment  is: 

<rank>  <n> 

<latl>  <lonl> 
<lat2>  <lon2> 


<latn>  <lonn> 

WDB-II  data  is  typically  read  from  tape.  Here  is  an  example  in 
which  the  Unix  command  dd  is  used  to  read  the  tape: 

dd  if-/dev/ntmtO  bs-=3600  conv=ascii  I  wdbtosff  >  eur.cil 

Written  by  Jim  Kilgore 

0.  S.  Naval  Research  Laboratory 

February  1989  */ 

♦include  <stdio.h> 

♦include  <sys/file.h> 

main (argc, argv) 
int  argc; 
char  **argv; 

( 

FILE  *of p=stdout , * f open ( )  ; 

register  int  linesegno, rank, if d-fileno (stdin) ,  npts, i; 

register  float  latdeg, latmin, latsec, longdeg, longmin, longsec, llatdeg, llongdeg 
char  ‘strncpy  O , ch2 ( 3) , ch3 (♦ ) , ns (2 ) , ew (2 ) , *pname-argv [  C ) ; 
char  ch20 (21),ch6(7),ch7(8),ch8(9); 

if  (argc  !«  1) 

{ 

(void)  fprintf  (stderr,  "usage:  %s:  <  infname  >  outfnameNn",  pname) ; 
return  -1;; 

) 

ch20[20]«' \0 ' ; 
ch6 [ 6) NO  '  ; 

Ch2 [ 2 !  —  * \ 0  '  ; 

ch7 [  7  ]  « ’ \0 ’ ; 

Ch8(8!='\0’; 
n  s  '  1  !  *  1  \  0  1  ; 
ew  J 1 ] = ’ \ 0  1  ; 

while  (1) 

( 

if  (read (lfd, ch20, 20)  —  0) 
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(voidj  fprintf  (stderr,  "%s:  EOF  while  reading  1st  ch20\n",  pname); 
return  -1; 


(void)  strncpy  ( 4ch7 [ 0 ] , ich20  !  0 ] , 7 ) ; 

Iinesegno=atoi (ch7) ; 

if  ((linesegno  !-  9)  Si  (linesegno  !  ■=  900000)  is  (linesegno  !=  999999)) 

( 

(void)  strncpy  (  Sch6 [ 0 j , Sch20 ( 9 ] , 6) ; 
npts=atoi (ch6) ; 

(void)  strncpy  (Sch2 (0 ] , Sch20 [7 ] , 2) ; 
rank=atoi (ch2) ; 

(void)  fprintf  (ofp,"%d  %dln", rank, npts) ; 

(void)  fflush  (ofp); 

for  (i=0;i<npts;i++) 

( 

if  (read (ifd, ch20, 20)  ==  0) 

{ 

(void)  fprintf  (stderr,  "%s:  EOF  while  reading  2nd  ch20\n",  pname) 
return  -1; 

) 

(void)  strncpy  ( Sch2 [ 0 ] , Sch20 ( 0 ] , 2 ) ; 
latdeg*=  (float)  atoi  (ch2)  ; 

(void)  strncpy  Uch2 (0i , Sch20 (21 ,2) ; 
latmin=(float) atoi (ch2) ; 

(void)  strncpy  ( 4ch2 (0 ) , Sch20 ( 4 ]  ,  2) ; 
latsec-(float) atoi (ch2) ; 

(void)  strncpy  ( ins [ 0 j , sch20 [ 6] , 1) ; 
llatdeg-latdeg  +  latmin/60 , +latsec/3600 . ; 
if  (ns  (0 ]  —  ’S') 

{ 

llatdeg=  -llatdeg; 

) 

(void)  strncpy  ( 4ch3 ( 0 ] , ich20 [ 1 ] , 3) ; 
longdeg- ( float ) atoi (ch3) ; 

(void)  strncpy  (4ch2[0),ich20[10],2); 
longmin- ( float) atoi (ch2 ) ; 

(void)  strncpy  (tch2 [0 ] , tch20 [ 12 ] ,  2) ; 
longsec- (float) atoi <ch2) ; 

(void)  strncpy  (iew (0] , i ch20 ( 14 | , 1 )  ; 
llongdeg=longdag+longmin/6G . +longsec/3600 .  ; 
if  (ew  (0  ]  —  ’W) 

( 

llongdeg-  -llongdeg; 

) 

(void)  fprintf  (ofp, "%f  %f\n", llatdeg, llongdeg)  ; 

(void)  fflush  (ofp); 

) 

) 

else 

( 

(void)  fprintf  (stderr,  "%s:  just  read  terminating  recordin’’,  pname); 
return  0; 

) 

) 

) 
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Appendix  B 


reduce  Program  Source  Code  Listing 

/*  REDUCE:  a  simple  data  reduction  algorithm  which  discards  points  i 
based  on  three  criteria:  minimum  distance  between  two  adjacent 
points,  maximum  distance  between  two  adjacent  points,  and  the 
angle  described  by  three  sequential  points. 

For  each  line  segment  (pi,  p2,  p3,  ...}  three  sequential  points  are 
considered:  pi,  p2,  and  p3 .  The  distances  between  pl-p2  and  p2-p3 
are  measured,  and  if  either  is  <»  a  minimum  distance  threshold 
then  point  p2  is  discarded.  Otherwise,  if  either  distance  is  >=  a 
maximum  distance  threshold  then  point  p2  is  discarded.  If  p2 
hasn't  been  discarded  yet  then  the  angle  A  described  by  the  three 
points  is  measured.  If  A  >-  max  angle  threshold  T,  where 
0  <-  T  <-  180  degrees  then  p2  is  discarded.  If  p2  was  discarded  in 
any  of  the  above  3  cases  then  p2  is  now  set  to  p3,  and  p3  is  set  to 
p4  .  Otherwise,  if  p2  was  retained  then  pi  is  set  to  p2,  p2  is  set 
to  p3,  and  p3  is  set  to  p4  .  The  process  begins  again  considering 
pi,  p2,  and  p3,  and  continues  in  this  manner  until  the  end  of  the 
segment  is  reached  for  all  segments. 

Written  by  Jim  Kilgore 
Naval  Research  Laboratory  */ 


#include  <stdio.h> 

♦include  <math.h> 

♦include  <sys/file.h> 

♦define  TRUE  1 

♦define  FALSE  0 

♦define  RADPERDEG  0.01745329 

♦define  PI  3.1415 

♦define  BUFLEN  30000 

♦define  disttxl,  yl,  x2,  y2)  (sqrt  ( (x2-xl) * (x2-xl) + (y2-yl) * (y2-yl) ) ) 
double  retainedlat (BUFLEN] , retainedlong (BUFLEN]  ; 


main (argc, argv) 
int  argc; 
char  **argv; 

{ 

FILE  *ifp=stdin, *ofp=stdout, *fopen () ; 

int  totalpoints, totalretained, discardp2, eof,rank,npts,i, nretained; 

int  mindistcrit,maxdistcrit, curvcrit; 

float  pllat , pi long, p21at , p21ong, p31at , p31ong; 

double  maxt hreshdeg , threshrad, at of  0 ,theta!23, sin()  ,  cos ()  ,sqrt ()  ,  acos ( )  ; 
double  minthreshdist,maxthreshdist; 
register  double  d, a, a2 , b, b2 , c, c2 , x; 
char  *pname=argv ( 0 ] ; 


< 


if  (argc 
< 

(void) 
ifname  > 
return 
1 


■  =  4) 

fprintf  (stderr,  "usage: 
of name\n",pname> ; 

-1; 


%s  maxthreshdeg  minthreshdist  maxthreshdxst 
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maxthreshdeg=atof (argv [  1  ] )  ; 

if  (maxthreshdeg  <  0.0  ! I  maxthreshdeg  >  180.0) 

( 

(void)  fprintf  (stderr,  "%s:  0.0  <=  maxthreshdeg  <-  180 . 0\n", pname) ; 
return  -1; 

) 

threshrad=maxthreshdec;*RADPERDEG; 

minthreshdist=atof (a-gv [2] )  ; 

if  (minthreshdist  <  0.0) 

( 

(void)  fprintf  (stderr,  "%s:  minthreshdist  must  be  >=  0 . C\n”, pname) ; 
return  -1; 

) 

maxthreshdist=atof (argv [ 3) ) ; 

eof-FALSE; 
totalpoints-0; 
totalretained«0 ; 

while  (!eof) 

{ 

if  (fscanf (ifp, "%d%d", trank, inpts)  ==  EOF) 

( 

eof-TRUE; 

) 

else 

{ 

(void)  fprintf  (stderr,  "%s:  rank=%d,  npts-%d\n", pname, rank , npts) ; 

totalpoints*totalpoints+npts; 

mir.distcrit-O; 

maxdistcrit-0 ; 

curvcrit-0; 

nretained=0; 

if  (npts  ««  1) 

{ 

/*  Read  point  but  reject  set:  */ 

if  ( f scanf (ifp, ”%f %f ", spllat , tpllong)  =  EOF) 

( 

(void)  fprintf  (stderr,  "%s:  unexpected  EOF  while  reading  pl\n", pname) 

return  -1; 

> 

) 

else  if  (npts  •-  2) 

( 

/*  Read  points  and  test  for  mindistcrit,  maxdistcrit:  */ 

if  (f scanf (ifp, "%f%f%f%f", Spllat, tpllong, tp21at, ip2iong)  ==  EOF) 

( 

(void)  fprintf  (stderr,  "%s:  unexpected  EOF  while  reading  pi,  p2\n",pn 

ame)  ; 

return  - 1 ; 

) 


16 


d=dist (pllong, pllat, p2 long, p21 at) ; 


if  (d  <=  minthreshdist ) 

mindistcrit=mindistcrit+2; 

else  if  (d  >=  maxthreshdist) 
maxdistcrit=maxdistcrit+2; 


else 

{ 

(void)  fprir.tf  (ofp,  "%d  %d\n",  rank,  npts)  ; 

(void)  fprintf  (ofp,  ''%f  »t\n%f  %f \n" , pllat , pllong, p21at , p21ong) ; 
nretained=2; 

) 


else 

( 

if  ( f scanf (ifp, "*f%f%f%f ", Spllat, spilong, 4p21at, ip21ong)  ==  EOF) 

( 

(void)  fprintf  (stderr,  "%s:  unexpected  EOF  while  reading  pi,  p2\n",pn 

ame)  ; 

return  -1; 

) 


/*  First  point  is  retained:  */ 

retainedlat [0] =pllat; 
retainedlong [0 ] =pllong; 
nretained++; 

for  (i-0;i<npts-2;i++) 

( 

if  (fscanf (ifp, "%f%f", 4p31at, 4p31ong)  ==  EOF) 

( 

(void)  fprintf  (stderr,  "%s:  unexpected  EOF  while  reading  p3\n",pnam 

return  -1; 

) 

a=dist (pllong, pllat, p31ong,p31at) ; 
b-dist (pllong, pllat, p21ong,  p21at)  ; 
c-dist (p21ong, p21at, p31ong, p31at)  ; 

if  (b  <=  minthreshdist  I  I  c  <=  minthreshdist ) 

( 

(void)  fprintf  (stderr,  "%s:  mindistcrit  violated\n", pname) ; 

discardp2“TRUE; 

mindistcrit++; 

) 


else  if  (b  >=  maxthreshdist  I  I  c  >=  maxthreshdist) 

( 

(void)  fprintf  (stderr,  "%s:  maxdistcrit  violatedXn" , pname) ; 

discardp2=TRUE; 

maxdistcrit++; 

) 

else 

f 

/*  Law  of  cosines:  a“2  *  b"2  +  c"2  -  2bc  cos  A:  */ 


a2=a*a; 
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b2=b*b; 

c2=c*c; 

x  =  (b2+c2-a2) / (2 . *b*c) ; 

/  Normalize  for  ocassional  roundoff  error  in  boundary  cases;  */ 

if  (x  <  -1.0)  /*  if  x  <  -1.0  it  is  just  SLIGHTLY  */ 

x=  -1.0;  /*  <  -1.0  and  should  be  corrected  to  -1.0  */ 

if  (x  >  1.0) 
x= 1,0; 

thetal23-acos (x)  ; 

/ *  Normalize  for  ocassional  roundoff  error  in  boundary  cases:  */ 

if  (thetal23  <  0.0)  /*  if  thetal23  <  0.0  it  is  just  SLIGHTLY  */ 

thet al23“0 . 0 ;  /*  <  0.0  and  should  should  be  corrected  to  0  */ 

if  (thetal23  >  PI) 
thetal23«PI; 

if  (thetal23  >-  threshrad) 

( 

discardp2=TRUE; 

curvcrit++; 

(void)  fprintf  (stderr,  1 %s:  curvcrit:  thenaI23=*f \n",  pname, theta 

123); 


else 

discardp2=FALSE; 

) 


if  (discardp2) 

( 

p21at=p31at;  /*  discard  p2  V 

p2 long-p31ong; 

) 

else 

( 

retainedlat [nretained] =p21at;  /*  retain  p2  */ 

retainedlong (nretained] =p2 long ; 

nretained++; 

pllat»p21at ; 

pllong“p21ong; 

p21at=p31at ; 

p21ong«p31ong; 

) 

) 


/*  Lastcheck:  Check  distance  between  last  retained  point  and  p3 :  */ 

p21at=retainedlat [ nretained-1 ] ; 
p21ong=retainedlong ( nretained-1 ] ; 

d=dist (p21ong, p21at, p31ong,p31at) ; 

if  (d  >»  maxthreshdist ) 

( 

(void)  fprintf  (stderr,  "%s:  lastcheck:  maxdistcrit  v i olation\n" , pname 

maxdistcrit++ ; 

) 


18 


( 

retainealat[nretained]=p31at; 
retainedlong[nretainedj  «p31ong; 
nretained++; 

i 


if  (.".retained  >2)  /*  Throw  out  sets  with  nretained=2  */ 

i 

/*  Write  rank,  nretained  to  output  file:  */ 

if  ( fprintf (ofp,  "%d  %d\n", rank, nretained)  ==  EOF) 

i 

(void)  fprintf  (stderr,  "%s:  fprintf:  could  not  write  rank,  nretained 
to  stdoutXn” , pname) ; 

return  -1; 

) 

for  ( i ■=  0 ;  i<nretained;  i++) 

( 

if  ( fprintf  (ofp,  *'%f  %f\n", retainedlat ( i ), ret ainedlong [ i ] )  ==  EOF) 

( 

(void)  fprintf  (stderr,  "%s:  fprintf:  could  not  write  retainedlat [ 
%d] ,  retainedlong ( %d]  to  stdout \n" , pname , i, i) ; 
return  -If 
) 

} 

} 

) 

totalretained-tocalretained+nretained; 

(void)  fprintf  (stderr,  "%s:  nretained-%d,  ", pname, nretained) ; 

(void)  fprintf  (stderr,  "pctretained=%f \n" ,( float ) nretained/ ( f loat) npts ) ; 
(void)  fprintf  (stderr,  "%s:  curvcrit=%d,  mindistcrit=%d,  maxdistcrit=%d\n 
\n" , pname, curvcrit , mindistcrit , maxdistcrit) ; 

) 

) 

(void)  fprintf  (stderr,  "%s:  done\n", pname) ; 

(void)  fprintf  (stderr,  "%s:  totalpoints-%d,  totalretained=%d,  pcttotalretai 
ned=%f \n", pname, total points, total  retained,  { float )totalretained/(fl oat )total point 
s)  ; 

(void)  fclose  (ifp) ; 

(void)  fclose  (ofp); 
return  0; 

) 
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Appendix  C 


bcor  Program  Source  Code  Listing 

/*  BCOR:  Corrects  boundary  crossing  errors 

Written  by  Jim  Kilgore 
Naval  Research  Laboratory  */ 


♦include  "bcor.h" 

♦define  MAXPCSNS  20000 

struct  position  p[MAXPOSNS]; 
extern  struct  _iobuf  _iob [ J ; 

main  large,  a rgv) 
int  arge; 
char  **argv; 

int  fsize,  eof  =  FALSE,  rank,  i,  npts,  cnt=0,  fstatO,  completed; 
long  ftell  () ,  fpos; 

FILE  *ifp=stdin, *ofp=stdout, *fopen () ; 
char  *pfname=argv[0) ; 
void  be ( ) ; 
struct  stat  buf; 

if  large  !=  1) 

( 

(void)  fprintf  (stderr,  "usage:  %s  <  ifname  >  ofname\n",  pfname) ; 
return  -1; 

) 

(void)  fstat  (fileno  (ifp),  tbuf ) ; 

fsize  =  buf.st_size; 

while  ( ! eof ) 

{ 

if  (fscanf (ifp, "%d%d", trank, tnpts)  »«  EOF) 

( 

eof “TRUE; 

) 

else 

( 

/* 

fpos  =  ftell  (ifp); 

completed  =  (int)  ((float)  fpos/  (float)  fsize  *  100.); 

(void)  fprintf  (stderr,  "\rcompleted :  td  pet",  completed); 

*/ 

if  (npts  >  MAXPOSNS) 

( 

(void)  fprintf  (stderr,  "\n%s:  td  >  MAXPOSNS:  MAXPOSNS  too  small\n" , pf na 
me ,  npt  s ) ; 

return  -1; 

) 

for  (i-0;i<-npts-l;i++) 
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if  (  f  scant"  (ifp,  "%f%f",  Sp[i]  .lat.  So  ;  i  ].  Ion)  ==  EOF) 

; 

(void)  fprintf  (stderr,  ”\n%s:  unexpected  EOF  while  reading  iat,  IcrAr. 

" , pf name) ; 

return  -I; 

) 

) 

(void)  be  (sent,  n.pts,  rank,  Sp[0),  cfp)  ; 

} 

i 

/* 

(void)  fprintf  (stderr,  "\rcorrections=%d  \n",  end; 

*/ 

return  0 ; 

) 


Program  Source  Code  Listing 

#include  "bcor.h" 


void  be  (ent,  n,  rank,  pos,  ofp) 
int  n,  rank,  *cnt; 
struct  position  pos  (  ) ; 

FILE  *ofp; 

( 

int  i,  j; 

register  float  d,  dxl,  dx2,  dx,  dy,  dyl; 

struct  position  pObound,  plbound; 

for  (i»l;  i<n;  i++) 

{ 

d  «  sqrt  ( (pos [ij ,lon-pos[i-l ) .Ion) * (pos [ i) .lon-pos [ i— 1 ] .ion) + (pos ! i ) .lat-po 
s  E 1-11 .iat) * (pos [i] . lat-pos ( i-1 ] . Iat) ) ; 

if  (d  >  MAXOIST) 

( 

/*  pObound  is  new  tail  of  first  segment,  plbound  is  head  */ 

/*  of  second  segment  */ 

if  (pos  !  i-1 )  . Ion  >  0  . ) 

< 

pObound. Ion  *  180.; 

plbound. ion  *  -180.- 


else 

( 

pObound. Ion  *  -180.; 
plbound. Ion  -  180.; 

) 


/*  Compute  boundary  point  by  simple  interpolation:  */ 

dxl  *  pObound. Ion  -  pos ( i-1 ). Ion; 
dx2  =  plbound. Ion  -  pos (i). Ion; 
dx  “dxl  -  dx2; 

dy  -  pos(i].lat  -  pos ( i-1 1 .lat; 
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dyl  =  cy/ax*dxl; 

pObound. iat  =  plbrund.iat  »  dyl  •*  pos  !  i-1 ;  .  Iat  ; 


/*  Write  rank,  n  points  in  first  segment:  */ 
(void)  fprintf  (ofp,  "%d  %d\n",  rank,  i  +  i); 


/*  Write  first  segment  to  file,  omitting  tail  point  pObound:  */ 
for  i j»0 ;  j< i ;  j  +  + ) 

(void)  fprintf  (ofp,  "%f  tf\r.",  pos! j]. iat,  pos  [  j  j  .  lor.)  ; 


/*  Now  write  tail  point  pObound  to  file:  */ 

(void)  fprintf  (ofp,  "%f  %f\n",  pObound. iat,  pObound. lor.)  ; 


/*  Set  pos! i-1 I  equal  to  head  point  plbound  of  second  segment,  */ 
/*  destructively  altering  old  value.  That's  ok  because  we  never  * 
/*  use  old  value  again:  */ 

pos ' i-1 5 , -iat  *  plbcund.iat; 
pos  i i-1 ). Ion  »  plbound. ion; 

**cnt;  /*  ent  i'  number  of  corrections  made  */ 


/*  Recursive  all  to  be  starting  with  pos[i-l]:  ’/ 
(void)  be  (ent,  n-i+1,  rank,  ipos[i-lj,  ofp); 

return ; 

) 


/•  Do  this  if  segment  passed  to  be  doesn't  need  to  be  split:  */ 

(void)  fprintf  (ofp,  "»d  %d\n”,  rank,  n) ; 
for  (i-C;  i<n;  i»+) 

(void)  fprintf  (ofp,  "»f  tf\n’’,  pos'ij.iat,  pos  I  i )  .  lor.)  ; 
return; 


bcor.h  Listing 


*  include 

*  include 

*  „  ucl ude 
« include 


<stdic.h> 
cm.ath .  h> 

<sys /types . h> 
<sy s /stat . h> 


•  define  TRUE!  1 

•  define  FALSE!  0 


♦define  MAXDIST  350. 


struct,  position 
( 
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flea'  lat; 
float  icn; 


Makefile  Listing 

CFLAGS  =  -f  s  itch 
LFLAG3  =  -im 


mai.i  :  bcor 

beer  :  bcor.c  bc.o 

re  S { CFLAGS }  bcor  .  bc.o 

bccr.c  :  bcor.c  bcor.h 

cc  ${ CFLAGS)  -c  bcor.c 

bc.o  :  bc.c  bcor.h 

cc  S (CFLAGS )  -c  bc.c 


-O  bcor  SiLFLAGS) 
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Appendix  D 


draw  Program  Source  Code  Listing 

/*  DRAW:  a  quick  ana  dirty  drawing  program 

cc  -f switch  draw.c  -c  dr  -L./lib  -lcore  -isunwindow  -lpixrect  -lm 

Note  that  lit-  and  include  directories  are  included  in  SPWD  because 
Su.nCcre  is  re  longer  supported  as  of  OS  4.C  release  and  those  files 
do  not  come  with  OS  distribution. 

Written  by  Jim  Kilgore 
Naval  Research  laboratory  */ 


♦include  "include /user core . h” 

♦include  <stdio.h> 

♦define  B1TLEN  lOCOOO 

int  pixwindd ( ) , bw2dd ( ) , cg2dd ( ) ; 

struct  vwsurf  vwsurf  *  DEFASLT_VWSL’RF  <bw2dd)  ; 

fioat  lat ; BUFLEN ) , ion [BUFLEN]  ; 

main  large, argv) 
int  arge; 
char  *’argv; 

( 

FILE  *if p, * f open  { )  ; 

float  aspe ct_ratio,r[125J,g[125)/b[12S), height, width ; 
double  iatmin,  lo.nmin,  latmax,  lonmax,  atof  O  ; 

^nt  eof , i , n, ir, ig, ib, rank, rankmax, npts, nf lies; 
char  ‘file,  *p,name=*argv  (  0]  ; 
void  cleanup  ( ) ; 

if  large  <  8) 

f 

(void)  fprintf  (stderr,  "usage:  %s  mIMIc  rankmax  iatmin  lonrin  latmax  lonmax 
fnamel  [fname2  f name3 . . . ’ \n", pname) ; 
return  -1; 

) 


rfiles  =  argc~"’ ; 

if  ( ’argv  '  1 j  ==  1 c 1 ) 

I 

vwsurf ,dd=cg2dd; 
height = 9CC ; 
w idt  h= 1 112; 


else  if  (*argv;ij  «*>  ’m1) 
l 

vwsurf ,dd=bw2dd; 
he iqht  =  900 ; 
width= 1 1 52 ; 

) 

else  if  (*argv[l)  —  ’MM 

( 
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vwsurf ,dd=bw2dd; 

height=1280; 

width=1600; 

j 

else 

< 

(void)  fprintf  (stderr,  "usage:  %s  blBlc  rankmax  latmin  lonm 
fnair.el  !fname2  f name3 . . . ] \n", pname) ; 
return  -1; 

) 

rankmax=*atoi (argv[2] ) ; 

latmin  «  atof  (argv[3)); 
lonmin  =  atof  <argv[4)); 
latmax  *  atof  (argv[5]); 
lonmax  -  atof  (argv(6)); 

/»  Fill  screen  for  quasi-eqruirectangular  projection:  */ 
aspect_r at io=height /width  ; 

/*  Gibberish  for  Core:  */ 

(void)  initialize_core (BASIC, NOINPOT,  TWOD)  ; 

(void)  initialize_view_surf ace (tvwsurf , FALSE)  ; 

(void)  select_view_sur f ace (tvwsurf ) ; 

(void)  set_ndc_space_2 (1.0, aspect_ratio)  ; 

(void)  set_viewport_2 (0.0, 1.0, 0.0, aspect_ratio) ; 

(void)  setwindow  (lonmin, lonmax, latmin, latmax) ; 

/*  Define  125  colors:  */ 

n = 0  ; 

for  (ir-0;ir<-4;ir++) 

( 

for  (ig-0; ig<-4 ; ig++) 

( 

for  ( ib=0; ib<=4 ; ib+*) 

( 

r [ n I = (float )  ir/4  . ; 
g I n ]  =  ( float) ig/4  .  ; 
r>(  n  i  =  (  float )  ib/4  .  ; 
n»+ ; 


) 

! 

(void)  def ine_color_indices (& vwsurf , 0, 124 , r , g, b) ; 

(void)  set _1 ' ne_index ( 12 4  )  ; 

(void)  set _ f i 1 i_ index ( 1 2 4 ) ; 

(void)  set  text  index  <  124)  ,- 

(void)  create_temporary_segment ( ) ; 

for  ( n  =  1 ;  n  <  =  r.  files  ;n  +  +) 

f ile=argv  f 6*n 1  ; 

if  (( if p*f open ( file, "r") )  NULL) 

( 

(void)  fprintf  (stderr,  "»s:  fopen:  could  not  open  \"%s\” 


in  latmax  lonmax 


for  read\n", pnam 
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e, file)  ; 

(void)  cleanupO; 
return  -1; 

) 

(void)  fprintf  (stderr,  "%s:  drawing  \"%s\”\.n", pname, file) ; 
eof =FALSE; 

while  (leaf) 

( 

if  (f scanf (ifp, "%d%d", Srank, Snpts)  ==  EOF) 
eof =TRUE; 
else 

{ 

if  ( npts  >  BUFLEN) 

( 

(void)  fprintf  (stderr,  "%s:  sorry,  npts-=%d  >  BUFLEN=%d\n", pname, npts 

, BUFLEN) ; 

(void)  cleanupO; 
return  -1; 

} 

for  (i=0; i<npts; i++) 

< 

if  ( f  scanf  ( ifp, "%f%f" , Slat [ i ], S Ion ( i ) )  ==EOF) 

t 

(void)  fprintf  (stderr,  "%s:  unexpected  EOF  while  reading  lat,  lon\n 

", pname) ; 

(void)  cleanupO; 
return  -i; 

>' 


if  (rank  <=  rankmax) 

( 

(void)  movc_abs_2 (Ion [0 ) , lat [ 0] )  ; 

(void)  polyline_abs_2 (Sion ( 1 ) , slat ( 1 ) ,  npts-1)  ; 
) 


) 

) 

(void)  f  close  <  ifp)  ,- 

) 

(void)  fprintf  ( stderr , "Hit  RETURN  to  clear  screen  and  guit : \n") ; 
(void)  getcharO; 

(void)  cleanup ( ) ; 
return  0; 

) 

void  cleanup)) 

( 

/*  Reset  colors  *  ' 

static  float  r[)  =  (1.0), 
g 1 ]  =  (1.0), 
b(!  =  (1.0); 

(void)  def ine_color_indices (Svwsurf , 0, 0, r,g, b) ; 

(void)  close_temporary_segment () ; 

/*  shut  down  core  V 

(void)  deselect_view_surface (svwsurf)  ; 

(void)  terminate_core ( )  ; 

) 
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Appendix  E 


suntokee  Program  Source  Code  Listing 

/*  SUNTOKEE:  converts  a  Sun  monochrome  (1  bit  deep)  standard 
rasterfile  (rasterfile  (5)  format  to  a  bitmap  file  in 
Xerox  Interlisp-D  format  which  is  compatible  with  KEE.  Sun 
color  rasterfiles  (6  bits  deep)  are  not  accepted  and  must  be 
first  converted  to  a  1-bit  deep  rasterfile  using 
rasf ilterStol (1) . 

Below  is  an  example  of  the  syntax.  Bitmapname  is  the  print 
name  that  will  be  given  to  the  bitmap. 

suntokee  bitmapname  <  rasterfname>  bitmapfname 

Written  by  Jim  Kilgore 

U.  S.  Naval  Research  Laboratory 

February,  1989 

Modified  8/9/90  to  correctly  read  and  convert  rasterfiles  whose 
ras_width  is  not  an  even  multiple  of  a.  */ 


tinclude  krasterf ile . h> 

♦include  <stdio.h> 

♦include  <sys/file.h> 

main  large,  argv) 
int  arge; 
char  **argv; 

{ 

FILE  *ofp=stdout, *fopen () ; 

int  lineno, pixelno, i,ifd-fileno(stdin) ; 

struct  rasterfile  rf; 

char  buf , buf hi, buf lo, *pname*argv [0 ) , *bitmapname=argv II]; 

if  large  !<=  2) 

I 

(void)  fprintf  (stderr,  "usage:  <  rasterfname  %s  bitmapname  >  bitmapf name\n 
,  pname) ; 

return  (-1 )  ; 

) 


/*  Read  the  header  and  print  it  out:  */ 

if  (read  (ifd, (char  *)  4 r f , sizeof ( r f ) )  ==  NULL) 

{ 

(void)  fprintf  (stderr,  "%s:  read:  unexpected  EOF  while  reading  header  stru 
turein",  pname) ; 
return  ( —  1 ) ; 

) 

/*  Magic  number  must  be  correct:  */ 

if  ( r f . ras_magic  !=  RAS_MAGIC) 

( 

(void)  fprintf  (stderr,  "%s:  wrong  magic  numberin'',  pname)  ; 
return  (-1 )  ; 

) 
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/*  ras_length  must  be  >  0:  */ 

if  ( rf . ras_length  <=  0) 

( 

(void)  fprintf  (stderr,  "%s:  bad  header  record,  rf . ras_length  must  be  >  C\n" 
,  pname)  ; 

return  (-1)  ; 


/*  Rasterfile  depth  must  be  1,  i.e.,  monochrome:  */ 

if  (rf ,ras_depth  !=  1) 

{ 

(void)  fprintf  (stderr,  "%s:  sorry,  rf.ras_depth  must  =  1 \n” , pname) ; 
return  (-1 )  ; 

) 

/*  Rasterfile  must  be  type  RT_STANDARD,  i.e.,  not  RT_BYTE_ENCODED, 
RT_FORMAT_RGB,  etc.  */ 

if  (rf  .ras_type  !**  RT_STANDARD) 

( 

(void)  fprintf  (stderr,  "%s:  sorry,  rasterfile  must  be  standard  format  image 
\n", pname) ; 

return  (-1); 

) 


lineno-0; 

pixelno=0; 

(void)  fprintf  (ofp, " (in-package  \'kee)\n"); 

(void)  fprintf  (ofp, " (def parameter  %s\n (expr-to-bitmap  \’(%d  %d  l\n\"", bitmapn 
ame,  r f . ras_width,  rf . ras_height) ; 


for  ( i=l ; i<=rf . ras_length; i++) 

( 

if  (readlifd, (char  *)  &buf,l)  ==  NULL) 

( 

(void)  fprintf  (stderr,  "\n%s:  unexpected  EOF  while  reading  image\n" , pname 

return  (-1 ) ; 

! 


/*  bufhi  =  upper  4  bits  of  buf;  buflo  =  lower  4  bits  of  buf:  */ 

bufhi  =  (buf  >>  4)  &  Oxf; 
buflo  =  buf  4  Oxf; 

/*  Translate  so  bufhi  and  buflo  lie  within  range  of  0x40  and  0x4f 
(between  ASCII  0  and  0) :  */ 


bufhi  +=  64; 
buflo  +=  64; 

(void)  fprintf 

pixelno  +=  8; 

if  (pixelno  >= 
ltiple  of  8  */ 

( 

*/ 

pixelno-0; 


(ofp,  ”»c%c", bufhi, buflo) ; 


rf.ras  width) 


/*  >=  because  rf . ras_width  may  not  be  even  mu 
/*  if  so,  screendump  pads  byte  with  zeros. 
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Iineno++; 

(void)  f print f  (stderr,  "\r%s:  completed:  %d  pet",  pname,  (int)  ((float)  1 
neno  /  (float)  r f . tas_height  *  100.)); 

(void)  fprintf  (ofp, "\""); 

if  (lineno  ==  rf . ras_height) 

(void)  fprintf  (ofp, ") ) ) \n") ; 

ej.se 

(void)  fprintf  (ofp, "\n\"") ; 

) 


(void)  close (ifd); 

(void)  fclose(ofp); 

(void)  fprintf  (stderr,  "\n%s:  done\n", pname) ; 
return  0; 
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Appendix  F 


Listing  of  File  box2 

1  5 

-90  -180 
90  -180 
90  180 
-90  180 
-90  -180 


Listing  of  File  box3 

1  5 
32  60 
23  60 
23  42 
32  42 
32  60 


Listing  of  File  box4 

1  5 

29.5  43 

30.5  47 

30.5  49 

29.5  49 

29.5  47 


